Can proper motions of dark-matter subhalos be detected? 
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One of the goals of NASA's Fermi Gamma-ray Space Telescope (formerly GLAST) will be de- 
tection of gamma rays from dark-matter annihilation in the Galactic halo. Theoretical arguments 
suggest that dark matter may be bound into subhalos with masses as small as 10 -4 -10 2 Af®. If 
so, it may be possible to detect individual subhalos as point sources in the Fermi Telescope. It 
has further been argued that some of these point sources may exhibit proper motions. Here we 
show that upper limits to the diffuse gamma-ray background constrain the number of subhalos close 
enough to exhibit proper motions to be less than one. 

PACS numbers: 95.85.Pw, 95.35.+d, 98.35.Gi 



I. INTRODUCTION 

One of the major goals of the recently launched Fermi 
Gamma-ray Space Telescope (formerly GLAST) [l| is to 
detect dark-matter annihilation into gamma rays, as ex- 
pected, with some flux, if dark matter is composed of 
weakly-interacting massive particles (WIMPs), such as 
the neutralino in supersymmetric models [2j. Theoret- 
ical arguments suggest that dark matter in the Galaxy 
may be distributed in subhalos [!, 0| with masses that 
may extend all the way down to 10 _4 -10 2 M e , if WIMPs 
make up dark matter 

References [(| pointed out that if WIMP substructure 
extends to such small scales, then it may be possible not 
only to detect individual subhalos as point sources with 
the Fermi Telescope, but also that up to ten of these 
point sources may exhibit proper motions. 

In this paper, we revisit this proposal. We use the 
diffuse gamma-ray background measured with the Ener- 
getic Gamma Ray Experimental Telescope (EGRET) Q 
to place an upper limit to the gamma-ray luminosity of 
these subhalos (see also Ref. [8]). This constraint then 
limits the number of detectable subhalos that are suffi- 
ciently close to exhibit proper motions. A similar argu- 
ment was also made in Ref. Q for a specific numerically- 
simulated subhalo model. Here, we present an analytic 
and fully general argument. The largest number of de- 
tectable proper motions is obtained if the subhalo mass 
distribution is sharply peaked around 1O _3 -1O _2 M0, 
and even then, the expected number of detectable proper 
motions is less than 1. 

Before proceeding with the detailed calculation, we be- 
gin in Sec. ITfl with an order-of-magnitude estimate. Sec- 
tion [TTT| then describes more precisely the EGRET upper 
limit to the diffuse background and the constraints im- 
plied for the subhalo population. We then evaluate in 
Sec. IIVI the number of detectable proper motions, first 
by assuming all subhalos have the same mass, and then 
by generalizing to the more realistic case of a power-law 
distribution of subhalo masses. In Sec. El we discuss the 
results and how they vary upon changing the assump- 



tions taken from their canonical values adopted in the 
main body of the paper. 



II. ORDER OF MAGNITUDE ESTIMATE 

If the Fermi Telescope can detect proper motions as 
small as 0.1° and subhalos have a transverse velocity 
^200 km s _1 , then a subhalo must be within a distance 
dpm ~ 0.35 pc to have a detectable proper motion over 
three years. If the point-source flux sensitivity is F sens , 
a source with a detectable proper motion must have a 
luminosity L ~ 47T(i 2 m i ;l scns . Each detected proper mo- 
tion thus implies a subhalo contribution to the luminosity 
density of order L/V ~ 3F SGns /G? pm . This implies a total 
luminosity, integrated over the volume of the Milky Way 
interior to the Solar radius, of Lmw ^ 47rrQF sens /d pm , 
where r Q = 8.5 kpc is our distance from the Galactic 
center (the ">" arises because the halo density most 
generally decreases with radius). Thus, a flux from the 
Milky Way of F MW ~ L M w/(47rr 2 ) ) > (r Q /d piri )F scns 
is implied for each proper motion detected. This is a 
large flux — larger than upper limits even for one detected 
proper motion — as the more detailed calculations below 
will show. 



III. EGRET CONSTRAINTS TO THE 
SUBHALO POPULATION 

Assume that a fraction / of the dark matter in the 
Galactic halo is composed of subhalos of mass M and 
luminosity L, 1 and we define a mass-to-gamma-ray- 
luminosity ratio T = M/L. Assume further that the 
radial distribution of subhalos follows the density profile 
p(r) of the Milky Way halo. 



1 Throughout this paper, the luminosity is the number (not en- 
ergy) of photons emitted per unit time; similarly, we deal with 
number fluxes and intensities. 
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The gamma-ray intensity I(ip) (units of pho- 
tons cm -2 s _1 sr -1 ), due to gamma rays from these sub- 
halos, is 



J(V) = 



4ttT 



dm P (r[liM) 



(1) 



along a line of sight at an angle ip from the Galactic 
center. Here, I is the distance along the line of sight ip; 



i.e., r 2 = + I 2 — 2r £cosV>- 



For p(r), we use the Navarro- Frenk- White (NFW) [l( 
profile, 



p(r) 



(r/r s )(l + r/r s ) 2 



(2) 



where p s = 5.4 x 10~ 3 Af Q pc -3 is the characteristic den- 
sity, and r s = 21.7 kpc is the scale radius. (These numer- 
ical values, taken to explain the virial mass and rotation 
curves [Tlj . are adopted for illustrative purposes, but the 
main conclusions do not depend on them.) This gives 
Pq = 7 x 1O~ 3 M pc~ 3 as the local density. 

The diffuse gamma-ray background measured with 
EGRET places an upper limit to /T _1 . If the subhalo 
gamma rays form a continuum spectrum, then we use the 
diffuse gamma-ray-background intensity measured with 
EGRET [l2j at ijj — 90° as a conservative upper limit 
to from subhalos (roughly the same bound obtains 
for ip = 180°). For line emission, on the other hand, we 
approximate the upper limit to the gamma-ray-line inten- 
sity, averaged over a 10° x 10° region around the Galactic 
center, by 2 x 10~ 6 (E/GeV)^ 1 ^ 2 cm~ 2 s _1 sr _1 over the 
energy range 0.1 GeV < E < lOGeV [13]. (Scaling our 
results to more conservative limits on line intensity [3| 
is trivial.) At an energy E — 10 GeV, the upper lim- 
its to both the continuum and line intensities turn out, 
coincidentally, to be 



/Tr- 



io 29 /, 



: Mq 1 B" 



(3) 



where !* < 1 is the fractional contribution from subha- 
los to the gamma-ray-background intensity, the total be- 
ing due, additionally, to traditional astrophysical sources 
and/or WIMP annihilation in the 1 — / fraction of the 
Milky Way halo that is smoothly distributed. We discuss 
further this assumption later on. 



IV. PROPER-MOTION DETECTION 



A. The mono-luminous case 

First, consider the case where the background is due to 
subhalos of a one mass scale M. We call this the mono- 
luminous case. The number of subhalos within a distance 
r of the Earth is 



iVsh(< r) 



47rr 3 fp 



3 M 
= 1.3 x 10 3 ^ 



o.i 



/ M 



\fM e 



(5) 



For the source to be detected at d pmi the luminosity 



l>m^sens; where i*sens 



L = Y -1 M has to be larger than 47rd 2 
is the point-source flux sensitivity; for the Fermi Tele 
scope, we take -F scns =10~ 10 !\o cm -2 s _1 at 10 GeV, 
with F w ~ 1 >5]. With the EGRET constraint to T 
[Eq. J5J], the condition for detection at d pm is M > M pm , 
where 



1.5 x 10- 2 fM F w I-%l 



(G) 



If M < Af pm , the subhalo can only be detected up to a 
distance [T- 1 Af/(4 7 rF scns )] 1 /2 < d pm . 

Thus, the number of detectable proper motions is 



N sh (< d pm ), 



N sh < 



T _1 M 



4irF aons 



1/2 



for M > M p 
for M < M, 



O.OgF^ 1 !,^^ x 



M 
M 



pmi 



for M > M n 



1/2 



, for M < M pn 



(7) 



To summarize, assuming the diffuse background den- 
sity measured by EGRET, and assuming that this back- 
ground is due to annihilations in subhalos, the maximum 
probability to detect proper motion in a nearby subhalo is 
^10%. This probability is maximized for subhalo masses 
M ~ 10 _2 /Mq. Larger-mass subhalos will provide a 
larger number of detectable point sources (that can be 
detected to larger distances), but fewer with proper mo- 
tions. On the other hand, there will be more nearby 
objects if they are less massive, but they will be too faint 
to be detectable. 



Suppose that a dark-matter subhalo moves transverse 
to the line of sight with a velocity 200 km s _1 . Then, in a 
three-year experiment, the subhalo will move an angular 
distance of 0.1° 6*o.i if it is at a distance 

d pm = 0.35pc^. (4) 

Here, we take the minimum proper motion detectable by 
the Fermi Telescope to be 0.1° (for 10 GeV gamma rays) 
as a canonical value [III, but keep the #o.i dependence 
in the following. Thus, a proper motion can be detected 
only if the subhalo is closer than d pm . 



B. Power-law mass and luminosity functions 

Now suppose, as before, that a fraction / of the mass of 
the halo is in subhalos, but now assume that the masses 
are distributed with a power-law mass function, 

* (r , M , . m^, m 

with power-law index —2, as suggested by numerical sim- 
ulations Here, A = M max /A! m i n , and Af m i n and 
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M max are the lower and upper mass limits. If the pop- 
ulation has a constant mass-to-luminosity ratio T, then 
the bound in Eq. ([3]) still applies. 

The differential number dN s h/dlnM of subhalos per 
logarithmic mass interval within a distance r from the 
Earth is given by Eq. ([5]), with the replacement / — > 
//In A. The simulations of Refs. [H suggest //In A = 
0.0064; typical values are In A = 35 and / « 0.2. The 
differential number of subhalos with detectable flux and 
proper motions per logarithmic mass interval is then 



dN, 



pm 



d\nM 



O.OO2F w 1 I*0q{ 



In A 
~35~ 



iE) > fOT M>M pu 



IT 1 ") V2 , f«»- ^1 - H 



(9) 



pm- 



Again the largest contribution comes from subhalos with 
masses around M pra . By integrating this expression over 
InM, we have 



pm 



(fe) ' 



for 



for M min < M pm < M max , 

(10) 



which is far smaller than in the mono-luminous case. 
This is easily understood: Detectable proper motions 
still come from subhalos with masses around M pm . How- 
ever, there are now fewer such subhalos, since the sub- 
halo masses are now distributed over a wider range. The 
contribution from M pm is thus smaller by a logarithmic 
factor In A. Although we focused only on the case of 
dn/dM cx M~ 2 here, the argument applies to any mass 
function. Again, the largest number of detectable proper 
motions comes when the mass distribution is sharply 
peaked around M pm . 



V. DISCUSSION AND CONCLUSIONS 

In this short paper, we revisit the detectability of 
proper motions of dark-matter subhalos. Our main re- 
sults are: (1) The EGRET diffuse-background constraint 
on the subhalo mass-to-luminosity ratio severely restricts 
the expected number of proper-motion detections with 
the Fermi Telescope to be less than 1. This argument 
does not rely on the details of any particle-physics model, 
in contrast to the earlier studies [f}]. (2) The number 
of detectable proper motions is maximized if the sub- 
halo mass distribution is sharply peaked around a mass 



M, 



pm 



0.01 /Mq. In this case, the probability to de- 
tect one subhalo with a proper motion is ~10%. (3) This 
maximal probability is obtained only when dark-matter 



annihilation from subhalos dominates the diffuse gamma- 
ray background (i.e., /* ~ 1), and it is reduced in propor- 
tion to the fraction 7* of the diffuse background due to 
subhalos. This fraction is likely to be much less than 1, as 
there may be significant contributions to EGRET's dif- 
fuse background from unresolved astrophysical sources, 
and perhaps from WIMP annihilation in a smoothly dis- 
tributed halo component. (4) If the subhalo mass func- 
tion follows a power-law mass function, the expected 
number of proper-motion detections is further suppressed 
by a logarithmic factor ln(M max /M m i n ), where M max and 
Mmin are masses of the largest and smallest subhalos. 
This factor could be ~30-40, in which case the chance 
of detection would be no greater than 1%. In this cal- 
culation, we have used a canonical gamma-ray energy of 
E = 10 GeV. The results can be scaled (but will not 
differ too much) by taking into account the E^ 1 / 2 de- 
pendence of 7* and a roughly similar energy dependence 
of the Fermi Telescope point-source sensitivity Fiq. 

Under these assumptions, the main result is that the 
probability of detecting proper motion in subhalos is less 
than one. It is natural then to ask under what conditions 
proper motions can be detected. There are two param- 
eters in our formalism which can be varied in order to 
explore this possibility. 

The quantity 7* can be used to parameterize the ef- 
fects of halo-modeling uncertainties. If p(r) increases 
less rapidly toward the Galactic center than in the NFW 
profile, then the predicted intensity will decrease and the 
bound in Eq. is weakened (for line radiation, which 
is bound by the Galactic center) accordingly. If, for ex- 
ample, we postulate an extreme scenario where the halo 
density p(r) is constant interior to the solar radius (hold- 
ing the local density fixed), then the bound in Eq. (T5|) is 
weakened by about a factor of 4. Likewise, 7* can also 
quantify the effect of a Galactocentric-radius-dependent 
fraction /. It may be that subhalos are more effectively 
disrupted closer to the Galactic center, and if so, then 
we should replace the fraction / in Eq. |T]) by a function 
f(r) that decreases monotonically as r — ► 0, and move 
f(r) into the integrand. If we guess f(r) = fo( r / r Q) a , 
we then find that for a = 1, the bound in Eq. ([3]) (with 
the replacement / — ► /o) is weakened by about a factor 
of 4 (the bound is weakened further by a factor of 6 if 
a = 2). Thus, if / has a radial dependence, then its ef- 
fects can be taken into account by allowing an increase in 
2* by a factor of 4 (for a = 1) to 6 (for a = 2). Note that 
this caveat applies only for line radiation (which we have 
bound by an upper limit to the gamma-ray flux from the 
Galactic center), but not to a continuum, for which we 
have used diffuse-background upper limits away from the 
Galactic center. 

The results are sensitive to the angular resolution and 
duration of the experiment. Throughout this calculation 
we assumed a threshold for proper-motion detection of 
0.1° and a timescale of three years. The number of de- 
tectable proper motions increases as 0q \ [see Eqs. 
and (|10p ]. Lowering the threshold for proper-motion de- 
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tection in a fixed timescale can be possible if a nearby 
subhalo is present, and detected at a high signal-to-noise 
ratio. In such a case the angular resolution of the instru- 

— 1/2 

ment will improve as ~ -ZV 7 , where iV 7 is the number 
of photons detected from the source. In addition, increas- 
ing the timescale of the experiment from three to, e.g., 
10 years (projected lifetime of the Fermi Telescope) will 
increase the probability of detection, as proper motion is 
linear with time. 

In summary, unless the distribution of subhalos within 
the solar radius is severely suppressed, and the Fermi 



Telescope operates for at least 10 years, the expected 
number of proper-motion detections is less than one. 
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